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Strong gravitational lensing: Why no central black holes? 
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Abstract. We investigate how central black holes (BHs) in galactic dark halos could affect strong gravitational 
lensing. The distribution of integral lensing probability with image separations are calculated for quasars of 
redshift 1.5 by foreground dark matter halos. The mass density of dark halos is taken to be the Navarro- Frcnk- 
White (NFW) profile such that, when the mass of a halo is less than 10 14 Mq, its central black holes or a bulge 
is included as a point mass. The relationship between the masses M, of supermassive black holes and the total 
gravitational mass Mdm of their host galaxy is adopted from the most recent literature. Only a flat ACDM model 
is considered here. It is shown that, while a single black hole for each galaxy contributes considerable but not 
sufficient lensing probabilities at small image separations compared with those without black holes, the bulges 
(which are about 100-1000 times larger in mass than a typical black hole) would definitely contribute enough 
probability at small image separations, although it gives too high probabilities at large separation angles compared 
with lensing observations. 
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1. Introduction 

Cold Dark Matter (CDM) has become the standard the- 
ory of cosmological structure formation. The ACDM vari- 
ant of CDM with fi m = 1 — Ha ~ 0.3 appears to be in 
good agreement with the available data on large scales 
(Primack |2002 ). On smaller (sub-galactic) scales, there 
seem to be various discrepancies, such: N-body CDM 
simulations which give cuspy halos with divergent pro 



files towards the center (Navarro, Frenk and White 1996 , 
|l997 , NFW hereafter); bar stability in high surface bright- 
ness spiral galaxies which also demands low-density cores; 
CDM models which yield an excess of small scale struc- 
tures; formation of disk galaxy angular momentum, which 
is much too small in galaxy simulations. Issues that have 
arisen on smaller scales have prompted people to propose 
a wide variety of alternatives to CDM, such as warm dark 
matter (WDM) and self-interacting dark matter (SIDM). 
Now that problems arise from galaxy-size halos and cen- 
ters of all dark matter halos, high-resolution simulations 
and observations are the final criterion. Recent highest- 
resolution simulations appear to be consistent with NFW 



Although further simulations and observations, inc ludin g 
measurement of CO rotation curves (Bolatto et al., 2002 ), 
may help to clarify the nature of the dark matter, it now 
appears that WDM and SIDM are both probably ruled 
out, while the small-scale predictions of ACDM may be in 
better agreement with the latest data than appeared to 
be the case as recently as a year ago. 

In addition to direct simulations and observations, 
gravitational lensing provides another powerful probe of 
mass distribution in the universe. Since mass within small 
scales only deflect light rays slightly, it is difficult to 
extract mass information from a single lensing event, 
and thus statistical gravitational lensing is needed even 
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Wambsganss, Cen, & Ostriker 1998; Porciani & Madau 
|2000| ; Keeton & Madau |2001[ ). Li & Ostriker $>002\ ) first 



(Klypin, 2002; Power et al. 2002) until scales smaller than 
about 1 kpc. Meanwhile, a large set of high-resolution op- 
tical rotation curves has recently been analyzed for low 
surface brightness (LSB) galaxies. One can also conclude 
that the NFW profile is a good fit down to about 1 kpc. 
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used the semi-analytical approach to analyze gravitational 
lensing of remote quasars by foreground dark halos in var- 
ious cold dark matter cosmologies. The mass function of 
dark halos they used is alternatively given by singular 
isothermal sphere (SIS), the NFW profile, or the gener- 
alized NFW profile. They found that none of these mod- 
els can completely explain the current observations: the 
SIS models predict too many large splitting lenses, while 
the NFW models predict too few small splitting lenses, 
so they proposed that there must be at least two pop- 
ulations of halos in the universe: small mass halos with 
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a steep inner density slope and large mass halos with a 
shallow inner density slope. The author conclude that a 
combination of SIS and NFW halos can reasonably re- 
produce the current observations. Similarly, Sarbu et al. 



(2001) investigated the statistics of gravitational lenses in 
flat, low-density cosmological models with different cosmic 
equations of state u>. It was found that COBE-normalized 
models with w > —0.4 produce too few arcsecond-scale 
lenses in comparison with the JVAS/CLASS radio survey, 
a result that is consistent with other observational con- 
straints on ui. 

When attention is attracted to alternatives of CDM 
dark matter density profile at small scales, another kind 
of dark matter — super-massive black holes in the cen- 
ters of most galactic halos is forgotten or ignored in this 
case, although the idea of detecting supermassive com- 
pact objects by their gravitational lensing effects was pro- 
posed very early (Press & Gunn 1973, Wilkinson et al. 
2001) and the lensing effects of Schwarzschild black holes 
in the strong field regime have been discussed in detail 
(e.g. Virbhadra & Ellis |2000| Frittelli, Kling & Newman 
200C; Bozza et al. 2001). On the other hand, cosmo- 



logical voids can form directly after the collapse of ex- 
tremely large wavelength perturbations into low-density 
black holes or cosmological black holes; such black holes 
can also be detected through their weak and strong lens- 
ing effects (Stornaiolo |2001| ). The observational evidence 
presented so far suggests the ubiquity of BHs in the nu- 
clei of all bright galaxies, regardless of their activity, and 
BH masses correlate with masses and luminosities of the 
host spheroids and, more ti ghtly, with stellar velocity 
dispe rsions (Magorrian e t al. 1998 ; Ferrarese & Merritt 
2000|; R avindran ath et al. |200l| ; Merri tt fc F errarese pOOla , 
200 lbt Wandel [2002| ; Sarzi et al. |2002| ). Most recent 
high-resolution observational data gives M./Mbuigc ~ 
10" 3 (Merritt & Ferrarese [2001cj ). Ferrarese ( |2002| ) fur- 
ther gave the relation between masses M. of supermassive 
black holes and the total gravitational mass of the dark 
matter halo in which they presumably formed 
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10 12 M, 







1.57 



(1) 



In this paper, we investigate the contributions of galactic 
central black holes to lensing probabilities at small im- 
age separations. Since ACDM cosmology and NFW profile 
are in good agreement with the available data of struc- 
ture formation on almost all scales as mentioned above, 
we only chose these two models respectively as cosmology 
and mass density function in our calculations. We model 
the lenses as a population of dark matter halos with an 
improved version of the Press-Schechter ( 1974 , PS) mass 
distribution function, and central BHs are considered for 
galaxy-size halos. 

The paper is organized as follows: the lensing equation 
is given in Sect. |2|, lensing probabilities arc calculated in 
Sect. 0, and discussion and conclusions are provided in 
Sect. |4[ 



2. Lensing equation 

The NFW profile is 



Pnfw — 



(2) 



r(r + r s ) 2 

where p s and r s are constants. We can define the mass of 
a halo to be the mass within r2oo (which is the radius of a 
sphere around a dark halo within which the average mass 
density is 200 times the critical mean mass density of the 
universe), 



Mi 



DM 



4tt 



pr 2 dr = 47rp s r 3 /(ci), 



(3) 



with ci = r2oo/r s the concentration parameter, the value 
of which is chosen to be 7/(l + z) (Bartelmann el al. 1998). 
And 



/(ci) 



ln(l + ci) 



ci 



(4) 



/ x(l + x) 2 x 1 + ci 

In flat ACDM cosmology, the constants p s and r s can then 



be expressed as (Li & Ostrikcr 2002), 



Ps = Peril [^m(l + zf + 
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3 /(d)' 
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(5) 



(6) 



where p cr u is the present value of the critical mass density 
of the universe, and Mi 5 is the reduced mass of a halo 
defined as M l5 = M DM /(lO 15 h- 1 M ). 

The surface mass density for NFW profile is 



J NFW 



(x) = 2p s r s x 



(l_ x 2)3/2 ! (A ^ / ! 



i(l/x) + Va; 2 -l-7r/2 
(x 2 -l) 3 / 2 



(x = l), 
(0 < x < 1 



Where x = \x\ and x = £/r s , £ is the position vector in the 
lens plane. The galactic central black holes arc assumed to 
be point masses, and we consider fist there is only a single 
black hole with mass M, for each galaxy. So the surface 
mass density for galactic halos each with a single central 
black hole can be written as 



galaxy 



(x) = M.S 2 {x) + E NFW (a;), 



(8) 



where 5 2 (x) is the two dimensional Dirac-delta function. 
The lensing equation with galactic central black holes con- 
sidered then is 



y = x- p s 



BH 



(9) 



where y = \y\, t] = yD^/D^ is the position vector in the 
source plance, in which D§ and are angular-diameter 
distances from the observer to the source and to the lens 
respectively. And 



Ps 



(10) 
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where E cr = (c 2 /AttG)(D^ /D£d£ s ) is the so called crit- 
ical surface mass density, in which D^ s is the angular- 
diameter distance from the lens to the source. And 



00) 



arctan \/ x 2 — 1 
Vx^l 



rctanh\/l — x 2 



(x > 1), 

(* = !). (11) 

(0 < x < 1). 



In Eq.(||), the term /bh stands for the contribution of a 
black hole, and by using Eq.(nl) and Eq.(0), it has the form 



/bh = 2.78 x 10- 4 /(ci)A/ 1 ° 5 



0.57 



(12) 



Since there is always more than one black holes in a bulge, 
and thus the bulge itself can act like a black hole, we can 
treat a bulge as a point mass in this paper as an extreme 
case. The mass of a black hole correlates linearly with 
that of its host bulge as M,/M hnlgQ ss 1CT 3 , so we can 
simply multiply the term / B h by 10 3 to stand for the 
contribution of a bulge. However, some light rays from 
the source will definitely travel across the bulge, so there 
must exist a kind of "effective" black hole with mass larger 
than a single "real" black hole but less than the bulge. 
In order to investigate the tendency of image separations 
contributed by different point mass, we can multiply /bh 
by, for example, 10 2 , etc. 

The lensing equations for three cases are plotted in 
Fig. [l] according to Eq. (|9]), where we have extended x 
and y to their opposite values because of symmetry. The 
full line, dashed line and dotted line, respectively, repre- 
sent the NFW lens with / BH = 0, 'NFW+BH' lens with 
/bh = 2.78 x 10- 4 /(ci)M0j 57 and 'NFW+bulge' with 
/bh = 2.78 x 10- 2 f(ci)M^ 57 , in which fi s = 0.49 and 
/(ci) = 0.91. 

We find that, as point masses, both a single central 
black hole and a bulge can more often produce small sep- 
aration images than the case when no central black holes 
are considered. This result will be further confirmed by 
the lensing probability given in next section. 

3. Lensing probability 

We choose the most generally accepted values of the pa- 
rameters for flat ACDM cosmology, for which, with usual 
symbols, the matter density parameter, vacuum energy 
density parameter and Hubble constant are respectively: 
n m = 0.3, fl A = 0.7, h = 0.75. 

The quasars of redshift z a — 1.5 are lensed by fore- 
ground CDM halos of galaxy clusters and galaxies, the 
lensing probability with image separations larger than A8 



is (Schneider, Ehlers, & Falco 1992) 



P(> A0) = 



dD L (z) 
dz 



dz 



n(M,z)a(M,z)dM. (13) 



Where D^(z) is proper distance from the observer to the 
lens located at redshift z 



Dl(z) = 



dz 



Ho Jo (l + z)y/n m (l + z) 3 + n A ' 



(14) 




0.2 



Fig. 1. Lensing equations: dashed line stands for a NFW 
lens with no central black holes; solid line stands for a lens 
with a NFW density profile plus a single central black hole; 
dotted line stands for a lens of NFW plus a bulge (treated 
as a point mass, an effective black hole). 



here c is the speed of light in vacuum and H is the current 
Hubble constant. The physical number density n(M, z) 
of virialized dark halos of masses between M and M + 
dM is related to the comoving number density n(M, z) by 
h(M, z) — n(M, z)(l + z) 3 ; the latter was originally given 
by Press & Schechter (1974), and the improved version is 



n(M, z)dM = j7f(M, z)dM, 



(15) 



where po is the current mean mass density of the universe, 
and 



f(M,z) 



2 6 c (z) din A 
7r MA din M 



exp 



5_M 
2A 2 



(16) 



is PS mass function. In Eq.(|l6|) above, A 2 (M) is the 
present variance of the fluctuations in a sphere contain- 
ing a mean mass M, 



1 r°° 

A 2 (M) = — / P(k)W 2 (k 



rM)k 2 dk, 



(17) 
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where P(k) is the power spectrum of density fluctuations, 
W{kryi) is the Fourier transformation of a top-hat window 
function 



W(kr M ) = 3 
and 



sin(fcrM) cos(fcrM) 
(fcr M ) 3 (kr M ) 2 



rM 



3M 

47T/90 



1/3 



(18) 



(19) 



In Eq.([l6|), 5 c (z) is the over density threshold for spherical 
collapse by redshift z (Navarro, Frenk, & White 1997): 



S c (z) 



1.68 



(20) 



where D{z) is the linear growth function of density per- 
turbation (Carroll, & Press |1992|) 



D(z) = 



g("(z)) 

g (n m )(i + zY 



in which 



, , 5 / 1 209x 



x 

140 



r 4/7 



and 



n m (i + z) ; 



i-n m + n m (i + z)3" 



(21) 



(22) 



(23) 



We use the fitting formulae for CD M power spectrum P(k) 
given by Eisenstein & Hu ( 1999 ) 

P(k) = AkT 2 {k), 



where A is the amplitude normalized to erg 
8/i- 1 Mpc) = 0.95, and 



with 

L = ln(e + 1.84g e ff), 



C = 14.4 



^n^Mpc- 1 ' 
325 



1 + 60.5^ 



l.n • 



(24) 
A(r M = 

(25) 

(26) 
(27) 

(28) 



We need to know the cross-sections in Eq.(^3[). Since 
we are interested in the lensing probabilities with image 
separations larger than a certain value AO (ranging from 
~ 10 arcseconds, for example), the cross-section is de- 
fined under the condition that the multiple images can 
be created. For the lenses with NFW profile, one can see 
from Fig^ that multiple images can be produced only 
if I V I < 2/cd where y cr is the maximum value of y when 
x < 0, which is determined by dy/dx = when /bh = 
in Eq.(||). For galaxy-size halos, the mass of which is con- 
fined to be less than 10 14 Mq through out this paper, a 



central black hole or a bulge as a point mass is consid- 
ered (see Eq.(||)). In this case, multiple images will always 
exist: when the source is close to the point caustic, i.e., 
when y is small, there are three images, two of which are 
within the Einstein circle and the third one the outside 
which Einstein circle; when y is large enough, there are 
two images, the weaker one is close to the center of the 
Einstein circle and the brighter one locates outside of the 
Einstein circle. So another condition is needed to define 
the cross-section, for which we use the brightness ratio 
between the brighter and weaker images just mentioned, 
and it is enough to set the ratio to be 10. 

The brightness ratio r for the two images is just the ra- 
tio of the corresponding abs olute values of magnifications 
(Schneider, Ehlers, & Falco |l992| ), 



M- 
where 

M+(y(z)) 



m-(s/(X)) = 



ydy 

x dx 

ydy 

x dx 



x>0 



x<0 



(29) 



(30) 



(31) 



Once the source position ?/ cr is determined by 

\fi+(Vcx)\ = 10|/i_(y cr )|, 



(32) 

the cross-section can be calculated, both with and without 
central black holes, as 



a(M, z) = iryl v r 2 s d{M - M mln ), 

where i)(x) is a step function, and M n 
lower limit of image separation 
r s Ax 2x r s 



A6 



(33) 

is determined by 
(34) 



and Eq.flq) as 



M min = 8.927xl0- 8 Af 15 (fi m (l + zf + fl A ) 



xq 



In Eq.(^J), we have approximated the image separation 
Ax to be 2xq, where Xg is the positive zero position of 
function y(x) , both when /bh = (for NFW lens only) 
and /bh (for galactic, NFW+BH/bulgc lenses) in 
Eq.(||), since image separation is insensitive to the source 
position y (Li & Ostriker, 2002 ). We plot the lensing prob- 
ability with image separations larger than AO in Fig. ^. In 
order to show the tendency of contributions to the lensing 
probability for different fraction of the bulge mass (the 
so called 'effective' black hole), we take the term /bh in 
Eq.(g) to be / BH = 2.78 x IQ- 1 f{ Cl )M^£ 7 , /bh = 2.78 x 
(5.0 x 10- 2 )/( Ci )M?b B7 and /bh = 2.78 x 10- 2 /(ci)M?b 6t ; 
they are represented by the first three lines from top down, 
respectively. When central black holes or bulges are in- 
cluded and treated as point masses, the mass of their host 
halos is confined to be less than 10 14 Mq, because Eq.(|lJ) 
strictly applies in the range 10 6 < M, < 2 x 10 9 Mq 
and 10 14 Mq is the upper-limit of galaxy mass (Ferrarese 
2002). 
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Fig. 2. Lensing probability with image separations larger 
than AO: the full line is for lenses of NFW density profile 
with no central black holes for halos at all scales, dot- 
ted line shows the case when a single central black hole 
for each galaxy-size halo is included. Other three lines 
show the cases when collectors of central black holes in the 
bulge are treated as an effective black hole. The dashed, 
dot-dash and dash-dot-dot-dot line, from top downwards, 
show respectively, the mass of the effective black hole at 
1000, 500 and 100 times that of a single 'real' black hole. 



4. Discussion and conclusions 

Our numerical results for lensing probability with image 
separations larger than A9 in five different cases are shown 
in Fig. |^. In all cases, lensing probabilities keep nearly con- 
stant until AO ~ 0.1 arc seconds, and obvious dropdown 
takes place at about 1 arc second if central black holes are 
included, which, of course, does not mean that the main 
lensing events have image separations larger than 0.1 arc 
seconds. As a matter of fact, in the NFW case (without 
galactic central black holes, the full line in Fig. ||), the 
lensing probability drops quite slowly in the whole range 
of image separations: AO ~ — 10 arc seconds; such a ten- 
dency would extend even to 30 arc seconds if it is plotted 
beyond this range, which implies a uniform distribution 
of lensing probability for its log value among image sep- 
arations. However, note that in the single black hole case 



(dotted line), the lensing probability drops to the same 
value of NFW at 2 arc seconds, which gives the influence 
range of a single black hole. In the range of ~ 0.1 arc sec- 
onds, the lensing probability for the single black hole case 
is about 3 times that for NFW. So, clearly, the contribu- 
tions from central black holes cannot be omitted, although 
such contributions alone are indeed not enough to explain 
the observational data. As we have pointed out, there is 
always more than one black hole in a galactic bulge, and 
the collector of black holes would make a bulge itself 'act 
like' a black hole. On the other hand, not all the mass 
of a bulge is concentrated in black holes, so if we treat a 
whole bulge as an extreme black hole, such a model would 
produce too many lenses at image separations larger than 
3 arc seconds compared with the JVAS / CLASS radio sur- 
vey. As mentioned above, we have sufficient reason to tune 
the fraction of a bulge mass to produce a 'right' profile of 
lensing probability at larger image separations required by 
observational data, but this 'sufficient reason' seems not 
make us produce sufficient lensing probabilities at smaller 
image separations, as shown by the dash-dot line in Fig. |[ 

However, we can attribute sufficient lens events at 
small image separations to galactic central black holes or 
the bulge. On the one hand, since this paper focuses on 
whether galactic central black holes would contribute con- 
siderably to the lensing probability, we have not consid- 
ered the effect of magnification bias, which would increase 
the final result provided here at all image separations. On 
the other hand, we have used an improved version of the 
PS halo mass function but not the 'best' version. The 
shape of the mass function predicted by standard PS the- 
ory (the improved version) is in reasonable agreement with 
what is measured in numerical simulations of hierarchical 
clustering from Gaussian initial conditions only for mas- 
sive halos,; less massive halos are more strongly clustered 
or less anti- biased than the standard PS predicted. Sheth 
& Tormen ( |l999| , ST) proposed a model that provides a 
reasonably good fit to the bias relation of less massive 
haloes as well as to that of massive halos. Note that cen- 
tral black holes are only found in galactic bulges,; ST's 
correction for mass function in the range of less massive 
halos would definitely change the lensing probability dis- 
cussed in this paper. Also note that one of the two im- 
ages produced by a galactic central black hole is close to 
the lens center and very faint; however, VLBI experiments 
can detect its existence (Hirabayashi 1990 ; Ulvestad 1999 ), 
and further radio lensing surveys would have the ability 
to identify high flux density ratio of the two images. How 
and to what extent lensing magnification bias, flux ratio 
and modified PS mass function may change the final result 
will be discussed in another paper. 
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